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Abstract 

In this paper, we consider the existence of a factorization of a monic, bounded motion 
polynomial. We prove existence of factorizations, possibly after multiplication with a real 
polynomial and provide algorithms for computing polynomial factor and factorizations. 
The first algorithm is conceptually simpler but may require a high degree of the polynomial 
factor. The second algorithm gives an optimal degree. 
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1. Introduction 

Let H[t] be the ring of univariate polynomials with quaternion coefficients, with the 
variable t commuting with the coefficients. The existence of factorizations of quaternion 
polynomials into linear factors is a classical result [2]. In [3], motion polynomials are 
defined as elements of DH[t] - the ring of univariate polynomials with dual quaternions 
- with real norms; these can be used to parametrize rational motions in Euclidean 3- 
space. The main result there is that a factorization into linear factors allows to construct 
a mechanical linkage that generates the desired motion. An adaption of the algorithm by 
[2] to the dual quaternion case indeed allows to factorize “generic” polynomials in IDHI [£], 
namely those whose primal part has no strictly real factors. For fixed degree, the set of 
generic motion polynomials is open and dense in the set of all motion polynomials. 

Since 2012, we have been wondering which non-generic motion polynomials do allow 
factorization into linear factors. One reason for our curiousity is a paradoxical fact: ratio¬ 
nal motions that are parametrized by generic motion polynomials have special properties, 
namely that their orbit curves have full cyclicity. The question is still not completely 
solved, but in this paper we give an affirmative answer for “bounded” motion polyno¬ 
mials. They always admit factorizations into products of linear rotation polynomials, 
possibly after multiplication with a real polynomial. This changes the motion polynomial 
but not the motion it parameterizes. Bounded motion polynomials are defined by the 
condition that the norm polynomials has no real roots. The kinematic meaning of this 
condition is that the orbits are bounded curves. It is also quite obvious that motions that 
can be generated by linkages with revolute joints (in particular, no translational joints) 
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are bounded, hence the results in this paper cover all cases for which there is a linkage 
with revolute joints. 

The results in this paper have been influenced by our study [6] of linkages producing 
straight line motions and by the factorization of planar bounded motion polynomials given 
in [!]■ 

The paper presents two different factorization algorithms for motion polynomials (Al¬ 
gorithms 3 and 4). Both compute a co-factor Q £ R[f] and a factorization of QM for 
a given bounded motion polynomial M £ BH[f]. The difference between these two algo¬ 
rithms is that Algorithm 3 is conceptually simpler, but it is calling Algorithm 2, introduced 
by [1] for the factorization of planar motion polynomials. Algorithm 4 is more complicated 
but it produces co-factors of optimal degree and does not depend on Algorithm 2. It is 
capable of producing non-planar factorizations of planar motion polynomials but can be 
specialized to yield planar factorizations as well. 

Structure of the paper. The remaining part of the paper is set up as follows. In Section 2, 
we recall the basic notations of dual quaternions and motion polynomials. Section 3 will 
focus on the factorization of a motion polynomial. We also give some concrete examples 
to support the algorithm. 

2. Dual quaternions and motion polynomials 

We start with a brief introduction to the dual quaternion model of rigid body dis¬ 
placements. In particular, we focus on one degree of freedom rational motions that can 
be parameterized by motion polynomials [7]. 

The dual quaternions form an eight-dimensional associative algebra over the real num¬ 
bers. This algebra is generated by the base elements 

1, i, j, k, e, ei, ej, ek 

and multiplication is defined via the relations 

i 2 = j 2 = k 2 = ijk = -1, e 2 = 0, ie = ei, je = ej, ke = ek. 

The set of dual quaternions is denoted by OH, the sub-algebra of quaternions H is gener¬ 
ated by 1, i, j, and k. Any dual quaternion can be written as h = p + eq with p, q £ H. 
The conjugate dual quaternion is h = p + eq and conjugation of quaternions is done by 
multiplying the coefficients of i, j, and k with —1. 

The norm of the dual quaternion h is defined as ||/i|| = hh. It equals pp + e(pq + qp) 
and is a dual number (an element of the sub-algebra D generated by 1 and e). A dual 
quaternion of norm 1 is called a unit dual quaternion. 

Dual quaternions have important applications in kinematics and mechanism science. 
This is due to an isomorphism between the factor group of unit dual quaternions mod¬ 
ulo ±1 and SE(3), the group of rigid body displacements. The rigid body displacement 
described by h = p + eq with hh = 1 maps the point x = x\\ + x- 2 j + a^k to 

pxp + pq — qp = pxp + 2 pq. 

Denote by MHI[£] the ring of polynomials in t with dual quaternion coefficients where 
multiplication is defined by the convention that the indeterminate t commutes with all 
coefficients. We follow the convention to write the coefficients to the left of the indetermi¬ 
nate t. Similarly, we denote by H[t] the sub-ring of polynomials with coefficients in H. The 
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conjugate polynomial to C = YH=o c *£i ^ MH[£] is C = Cjt* and the norm polynomial 
is CC. Its coefficients are dual numbers. If C = c *£% the value C{h) of C at h £ OH 

is defined as C{h) = 5Z"=o We a ^ so define C( oo) := c n . 

A polynomial M = P + eQ £ KMHI[£] is called a motion polynomial if PQ + QP = 0 and 
its leading coefficient is invertible. Usually we will even assume that the leading coefficient 
is one (the polynomial is monic). This can be accomplished by left-multiplying M with the 
inverse of the leading coefficient and often constitutes no loss of generality. The defining 
conditions of a motion polynomial ensure that its norm polynomial has real coefficients. 

A motion polynomial M = P + eQ acts on a point x = aqi + x%j + 2 : 3 k according to 


x i->- 


PxP + 2 PQ 


( 1 ) 


This equation defines a rigid body displacement for all values t 6 8 U { 00 } that are not 
zeros of P. Any map of the shape (1) with a motion polynomial M = P + eQ is called 
a rational motion. We also say that the motion polynomial parameterizes the rational 
motion. The motion’s trajectories (orbits of points for t £ R U { 00 }) are rational curves. 
It is known that any motion with only rational trajectories is parameterized by a motion 
polynomial [5]. 

The simplest motion polynomials are of degree one and can be written as M = t — h 
where h. — h £ 1 and hh £ R. They parameterize either rotations about a fixed axis or 
translations in a fixed direction. We speak of the rotation or translation quaternion h and 
the rotation or translation polynomial t — h, respectively. In this paper we are concerned 
with the factorization of motion polynomials into the product of rotation polynomials. 
These are distinguished from translation polynomials by having a primal part not in R[f]. 


3. Factorization 

In [3] it has been shown that a generic monic motion polynomial M = P+eD of degree 
n admits factorizations of the shape 

M = - h n ) (2) 

with rotation polynomials t — hi ,..., t — h n . Here, the term “generic” means that the 
primal part P of M has no real factors. The factorization (2) can be computed by the 
non-deterministic Algorithm 1. The details of this algorithm are explained in [3] but some 
comments are appropriate at this place. 

• In all our algorithms, we denote concatenation of lists by the operator symbol “+”. 
List concatenation is not commutative: The list L\ + L 2 starts with the elements of 
L 1 and ends with the elements of £ 2 - 

• By genericity of M, the norm polynomial PP is real and positive. Hence, it is the 
product of n quadratic, real factors which are irreducible over R. 

• The choice of a quadratic factor in Line 5 is arbitrary. Different choices result in 
different factorizations. In general, there are n! factorizations of the shape (2), each 
corresponding to a permutation of the quadratic factors of PP. 

• For left polynomials with dual quaternion coefficients in our sense, right division is 
possible: Given two polynomials M, N £ H[£] with N monic, there exist unique 
polynomials Q, R £ H[f] with M = QN + R and deg I? < deg N. 
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• The dual quaternion hi in Line 6 can be computed as zero of the linear polynomial 
Ri obtained by writing M = QMi + Ri (polynomial division). The assumptions on 
M guarantee existence of a unique zero over the dual quaternions but the algorithm 
may fail at this point if these assumptions are not met. 

• We may exit the algorithm after just one iteration to find a linear right factor of M, 
that is, write the motion polynomial as M = M'(t — h). This we will often do in our 
factorization algorithm for non-generic motion polynomials. 


Algorithm 1 GFactor 

Input: M = P + eD G OH [t]. a monic, generic motion polynomial of degree n. 
Output: A list L = [Li ,..., L n \ such that M = L\ - ■ ■ L n . 


1: L <- [] 

2: F 4— [Mi, . M n \ 

3: 

4: repeat 

5: Choose Mi G F and set F <— F 

6 : Compute hi such that = 

7: L 4 — [t — hf\ T L 

8 : M 4— M/{t — hi) 

9: until degM = 0. 

10 : Return L = [Li, L 2 ,..., L n ]. 


> (empty list) 
> Each Mi G M[i], i = 1,..., n is a 
quadratic, irreducible factor of PP G R[t]. 


- [Mi]. 

M{hi) = 0. 

> (add t — hi to start of list) 
> (polynomial division) 


For later reference, we state the result of [3, Theorem 3] as a lemma. We do this in 
a form that highlights the dependence of the factorization on an ordering of the norm 
polynomial’s quadratic factors. 

Lemma 1. Given a generic, monic motion polynomial M of degree n with MAi = 
Mi ■ ■ ■ M n and monic, quadratic and irreducible real polynomials Adi, ■ • •, M n , there exist 
rotation quaternions hi ,..., h n such that M = (t — hi) ■ ■ ■ (t — h n ) and Mi = (t — hi){t — hi) 
for i = 1,... ,n. Different labeling of the quadratic factors of MAI give different factor¬ 
izations. 

Here are examples of non-generic motion polynomials with exceptional factorizations. 

Example 1. The motion polynomial M := t 2 + 1 + ei is not generic. A straightforward 
computation shows that no linear motion polynomials t — hi and t — h -2 in OH[t] with 
M = (t— hi){t — / 12 ) exist. The motion parameterized by M is a translation with constant 
direction. 

Example 2. Non-generic motion polynomials with infinitely many factorizations exist. One 
example is M := t 2 + 1 — et(it — j). It can be factorized as M = (t — hi)(t — hf) where 

hi = k — e(ai + (6 - l)j), h 2 = -k + e(ai + fej) 

and a, b are arbitrary real numbers. The motion parameterized by M is a circular trans¬ 
lation. Any of the infinitely many factorizations of M corresponds two one leg of a paral¬ 
lelogram linkage that can generate this motion. 
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Example 3. The motion polynomial M := t 2 — (1 +j)f +j — e((i + k)£ — 2k) can be factored 
as 

M = (t — 1 — ei)(t — j ek) = (t — j — e(i + 2k ))(£ — 1 + ek). 

The polynomial factors t — 1 — ei and t — 1 + ek parameterize, however, translations, 
not rotations. The reason for this is the possibility to factor the primal part of M as 
t 2 — (1 + j )t + j = (t — l)(t — j). For t = 1, the motion parameterization becomes singular 
and the trajectories pass through infinite points. 

We will present a method to factor even the motion polynomials of these examples into 
products of linear rotation polynomials. This will be made possible by allowing alterations 
of the given motion polynomial that change its kinematic and algebraic properties in an 
“admissible” way. This alterations are: 

1. Multiplication of M with a strictly positive real polynomial Q and factorization of 
QM instead of M. This is an admissible change because M and QM parameterize 
the same motion. This “multiplication trick” has already been used in [1] for the 
factorization of planar motion polynomials. 

2. Substitution of a rational expression R/Q with R, Q £ R[t] for the indeterminate 
t in M and factorization of Q desM M(R/Q) instead of M. This amounts to a not 
necessarily invertible re-parameterization of the motion. In particular, it is possible 
to parameterize only one part of the original motion. 

Multiplication with real polynomials does not change kinematic properties but gives 
additional flexibility to find factorizations in otherwise unfactorizable cases. In order to 
explain the meaning and necessity of substitution of real polynomials, we first give an 
important definition. 

Definition 1. A motion polynomial M = P + eD is called bounded, if its primal part P 
has no real zeros. 

Generic motion polynomials are bounded. Bounded motion polynomials parameterize 
precisely the rational motions with only bounded trajectories. If the motion polynomial 
is not bounded, zeros of the primal part belong to infinite points on the trajectories. 
For this reason, unbounded motion polynomials can never be written as the product of 
linear rotation polynomials. For example, we can never succeed in finding a factorization 
(t — h\)(t — h, 2 ) with rotation quaternions h±, /12 of the motion polynomial in Example 3 
as it has unbounded trajectories. 

Unbounded motion polynomials can always be turned into bounded ones by an ap¬ 
propriate substitution. This is the reason, why we henceforth restrict our attention to 
bounded motion polynomials. The kinematic meaning is that only a certain portion of 
the original trajectories is actually reached during the motion. Finally, we assume that 
our motion polynomials are rnonic. This is no loss of generality. If M is bounded, the 
leading coefficient c n of M is invertible and we may factor c~ x M instead. This amounts 
to an admissible change of coordinates. 

To summarize and give a precise problem statement: Given a bounded, rnonic motion 
polynomial M, we want to find a real polynomial Q and a list of linear rotation polynomials 
L = [t — hi ,..., t — h n \ such that QM = (t — hi) •■•(£ — h n ). In this case we say that “M 
admits a factorization”. We will not only prove existence of Q and L , we will also provide 
a simple algorithm for computing appropriate Q and L , provide a bound on the degree 
of Q (and hence also on the number of polynomials in L) and present a more elaborate 
algorithm that produces a polynomial Q of minimal degree. 
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3.1. Factorization of non-generic cases 

On particular case for which existence of factorizations of non-generic motion polyno¬ 
mials has already been proved to exist is planar kinematics [ 1 ]. 

Definition 2. A motion polynomial M is called planar, if it parameterizes a planar 
motion (a subgroup consisting of all rotations around axes parallel to a fixed direction and 
translations orthogonal to that direction). 

Examples of planar motion polynomials are obtained by picking coefficients in (1, i, ej, ek). 
In [1], the authors showed that for every monic, bounded, planar motion polynomial M 
of degree n a real polynomial Q of degree deg Q < n exists such that QM admits a fac¬ 
torization of the shape (2). Input and output of this planar factorization algorithm are 
displayed in Algorithm 2. We list this algorithm only for the purpose of later reference. 
For details we refer to [1]. 


Algorithm 2 PFactor (planar factorization algorithm of [1]) 

Input: M = P + eD G a planar, bounded, monic motion polynomial. 

Output: Q G R[t], list L = [L\, L 2 , ■ ■ ., L n \ of linear rotation polynomials such that 
QM = L 1 L 2 ■ • ■ L n . 


The first factorization procedure we propose is of theoretical interest. It is displayed in 
Algorithm 3. It is based on the algorithm for factorization of planar motion polynomials 
and produces a real polynomial Q and a factorization of QM for a monic and bounded but 
not necessarily generic motion polynomial M. It is conceptually simpler than Algorithm 4 
below but non optimal as far as minimality of deg Q is concerned. In its listing, we denote 
by GRPF(M) the greatest real polynomial factor of a quaternion polynomial M G H[t]. 
Lines 2 to 5 of Algorithm 3 are based on the factorization 

MTT = (R\T + eD)TT = (R^T + eDT)T 

of MTT into the product of a planar motion polynomial and a polynomial T G H[f]. 
Together with [1], Algorithm 3 proves existence of a factorization: 

Theorem 1. Given a bounded, monic motion polynomial M G BEI[t] there always exists a 
real polynomial Q such that QM can be written as a product of linear rotation polynomials. 

3.2. Factorizations of minimal degree 

Now we should further elaborate on the minimal possible degree of the real factor 
Q that makes factorization possible. In the planar case, Algorithm 2 gives the bound 
degQ < deg M and this bound is known to be optimal [ 1 ], The upper bound achievable 
with Algorithm 3 is worse. Let m = deg Ad and r = deg R\. Then, the degree of <5 in Line 4 
is bounded by 2(m — r) and the degree of P in Line 4 is bounded by r + 2(m — r) = 2m — r. 
Hence, the degree of Q at the end of Algorithm 3 is bounded by 2 (to — r) + 2(2 m — r) = 
6 to — 4r. Because of r > 2, this gives the bound degQ < 6 m — 8. However, also in the 
spatial case the bound deg Q < deg M holds true. This is guaranteed by Algorithm 4. 
Here are a few remarks on Algorithm 4. 

• In Algorithm 4, we mainly treat the case where the primal part P of the motion 
polynomial M = P + eD has a non-constant real factor f?i = GRPF(P). Otherwise, 
we just resort to factorization of generic motion polynomials (Algorithm 1 ). 
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Algorithm 3 Factorl 

Input: M = P + eD £ BH[f], a monic, bounded motion polynomial with real quadratic 


factor in its primal part, Q £ M[f], list L of linear motion polynomials. Initially, Q = 1 


and L = [] (empty list). 


Output: Q and L = [Pi, L 2 , ■ ■ ■, L n ] such that QM = P 1 P 2 • • • L n . 

1 

Write P = RiT where Pi = GRPF(P). 


2 

If deg T / 0 Then 


3 

L 4 — L + GFactor(T) 

> Append linear factors of T to P. 

4 

Q •<— TT, P t— R\TT , D •<— DT , and M<-P + eD 

5 

End If 


6 

Factor MP = (P + e(P>ii + £> 2 .] + D 3 k))P = 

(P + eP>ii)(P + eP> 2 j + eD 3 k). 

7 

Qi, Pi = PFactor(P + ePhi) 

> (planar factorization) 

8 

Q 2 ,L 2 = PFactor (P + eP) 2 j + eD 3 k) 

> (planar factorization) 

9 

Q QQ 1 Q 2 = QP 2 

> (because Qi = Q 2 = P) 

10 

L 4 — L 2 + L\ + L 

> Concatenate lists of linear factors. 

11 

Return Q , P 



• The complexity of a monic bounded motion polynomial M = P + eD £ OH[f] in 
Algorithm 4 is a triple of integers 

comp(M) := (a, £, 7 ), 

a := deg(gcd(P, P, DD)), 

/3 := deg(gcd(P,P)), 

7 := deg(P), 

where deg(a) is the degree of the polynomial a and gcd(a, b) £ K[t] is the greatest 
real common factor of polynomials a and b. With this definition, gcd(a,a) is the 
greatest real polynomial factor of a. In each step of the recursive Algorithm 4, we 
try to construct M' such that comp(M') < comp(Af) with lexicographic order, e.g., 
(4, 2, 5) < (4,4, 3), (4, 2, 2) < (4, 2,3). Then we recursively call FactorAll with M' 
as argument. As soon as a = /3 = 7 = 0, Algorithm 4 terminates. 

• The computation of quaternions hi and h r in Lines 33-34 is based on Lemma 1 and 
[4, Theorem 3.2]. One of this theorem’s statements is that the set of quaternion 
roots of the irreducible quadratic polynomial Q = t 2 + bt + c £ R[f] is 

{^( - b + \/4c - b 2 (xii + X2$ + 2:3k)) | (xi,xi,x 3 ) £ S' 2 } (3) 

where S 2 is the unit 2-spliere in R 3 . In particular, for every unit vector (x 3 ,X 2 ,x 3 ) £ 
S 2 , there is a the quaternion root q whose vector part is proportional to x\i + X 2 J + 
x 3 k. Also note that Q = (t — h)(t — h) if h is a quaternion root of Q. In the 
algorithm, we can pick an arbitrary zero h r of Pi and compute D r by polynomial 
division. Then we compute hi as zero of the remainder polynomial R in the division 
D = QM + R with M = (t — h r )(t — h r ), as in one iteration of Algorithm 1 and Di 
again by polynomial division. 

• The computation of quaternions hi and h r in Lines 10-12 and Lines 24-26 of Al¬ 
gorithm 4 is again based on Lemma 1 but also on Lemma 2 below. Consider, for 
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Algorithm 4 FactorAll 

Input: M = P+eD £ OH[f], a monic, bounded motion polynomial of complexity (a, /3, 7 ), 
Q £ R[i], lists Li , L r of linear motion polynomials. Initially, Q = 1, Li = [], L r = []. 
Output: Q , Li , such that with Li + L r = [Li, L 2 , ..., L„] we have QM = P 1 L 2 ■ ■ • L n . 

1 : If P has no real factors Then 
2: Return Q, L/, L r + GFactor(M). 

3: End If 

4 : Let R\ be the GRPF of P, i.e., P = PiT, degP = f3. 

5: Let a := deg(gcd(P, P, DD)) = deg(gcd(Pi, DD)). > comp(M) = (a,/3, 7 ). 

6 : If gcd(Pi,PP) = 1 (a = 0) Then 
7: If gcd(Pi,PT) = 1 Then 

8 : If T = 1, i.e., P is real Then 

9: Let Pi be a quadratic real divisor of P, i.e., P = P\P'. 

10 : Compute quaternion roots h r , hi of Pi such that 

11 : hi h r , D(t — h r ) = (t — hi)D', > (Lemma 1, Lemma 2) 

12: (t — hi)D'{t — h r ) = DP\. 

13: Q 4— QP\, Li 4 — Li + [f — ftj], L r <— [i — /i r ] + L r , 

14: M't— P'(t — hi)(t — h r ) + eD' . t> comp(M') = (0, — 2, 7 ). 

15: Return FactorAll(Af', Q, Lp L r ) 

16: Else 

17: Let Pi be a quadratic real divisor of TT. 

18: Compute a common zero h of Pi and M such that 

19: Pi = (t — h)(t — h), M = M'(t — h). > comp(AL') = (0, /3 ,7 — 1). 

20: Return FactorAll(M', Q, Li, [t — h\ + L r ) 

21 : End If 

22 : Else 

23: Let Pi be a quadratic real divisor of gcd(Pi,TT), i.e., P = P'P\. 

24: Compute quaternions roots h r , hi of Pi such that 

25: Pi{h r ) = 0, T(hr) 7 ^ 0, T(hi) 7 ^ 0, > (Lemma 1, Lemma 2) 

26: DP\ = D(t — h r )(t — h r ) = (t — hi)D'{t — h r ). > (Lemma 1, Lemma 2) 

27: <3 4— QPi, Li 4— L; + [t — h{\ , L r 4— [t — /l r ] + L r , 

28: M' 4— (t — hi)P'(t — h r ) + eD'. > comp(M') = (0,/3 — 2 , 7 ). 

29: Return FactorAll(M', Q, L;, L r ) 

30: End If 

31: Else (a > 2) 

32: Let Pi be a quadratic real divisor of gcd(Ri, DD). 

33: Compute quaternion roots h r , hi of Pl such that c> (Lemma 1) 

34 : D = (t — hi)Di = D r (t — h r ) and P = (t — hi)Pi = P r (t — h r ) > (Lemma 1) 

35: If deg GRPF(P) < deg GRPF(P r ) Then 

36: Li 4— L; + [t — /q] , M 7 4— P; + eD;. > comp(M') = (a — 2, /? — 2 ,7 — 1). 

37: Return FactorAll(M', Q, Li, L r ) 

38: Else 

39: L r 4— [t — h r ] + L r M' i — P r + eD r . > comp(M') = (a — 2 , /?, 7 — 1 ) 

40: or comp(M') = (a — 2, (3 — 2 ,7 — 1). 

41: Return FactorAll(M', Q, Li, L r ) 

42: End If 

43: End If 






example, the situation in Lines 10-12. We may prescribe h r arbitrarily as a root of 
Pi, see (3). Then we use polynomial division (over ME) to find Q,R £ H[f] such 
that ( t — h r )D = QPi + R and compute hi as unique zero of the linear remainder 
polynomial R. Using polynomial division once more, we then find D' such that 
(■ t-h r )D = D 7 (t-h ~ l ). 

Lemma 2 . Let Q 6 R[l] be a quadratic polynomial that is irreducible over M, D £ i[(] a 
polynomial with gcd(DD,Q) = 1 and O the set of quaternion roots of Q. Then the map 
fQ,D' O —> O, hi i-> h r with h r being the common root of (t — hi)D and Q is a well-defined 
bijection. Moreover, fQ,D{h) ^ h for all h £ O. 

Proof. Our proof is based on results of [3] that state that the quaternion roots of a poly¬ 
nomial P £ H[t] are also roots of the quadratic factors of PP. Moreover, h is a root of P 
if and only if t — h is a right factor of P [3, Lemma 2]. 

By (3), the set O is not empty. The norm polynomial of (t — hi)D has the quadratic 
factor Q. Hence, there exists a quaternion root h r £ O of (t — hi)D. This root is unique 
because of gcd (DD, Q) = 1 and the map /q,d is well-defined. 

If fQ,D (h) = h for some h £ O, there exists D' £ H[t] with D = (t — h)D'(t — h) and 
we get a contradiction to gcd {DD, Q) = 1: 

DD = (t — h)D\t — h)(t — h)ly(t -h) = Q(t - h)D'W(t - h ). 

By a linear parameter transformation 1 1 ->- at + b with a, b £ R we can always achieve 
that Q is a real multiple of t 2 + 1. Hence, it is no loss of generality to assume Q = t 2 + 1 
when proving bijectivity of /q,d. Using polynomial division we find K £ H[f] and a, b £ R 


with D = K(t 2 + 1) + at + b. Then we have 

(t - hi)D = (t - hi)K(t 2 + 1) + (t - h){at + b) 

= ((£ — hi)K + a)(t 2 + 1) + (b — hia)t — a — hib. 

As already argued, there is h r = fo(hi) £ O such that 

(b — hia)h r — a — hib = 0. (4) 

If there is h[ ^ hi with /d(/i() = h r then we also have 

(b — h'^hr — a — h[b = 0. (5) 

Subtracting Equations 4 and 5 yields 

{h[ — h{)ah r + (h[ — hi)b = 0 . ( 6 ) 

As h[ — hi ^ 0, we have ah r + 6 = 0 and this implies D = K(t 2 + 1) + aft — h r ). But 
then deggcd(PP, t 2 + 1) >0 would contradict our assumptions. Hence fo is injective. To 
prove surjectivity, observe that for any h r £ O , there is hi such that (t — h r )D = D'(t — h{) 
by injectivity of fjj. But then we have foihi) = hr. □ 

The termination of Algorithm 4 is guaranteed by the following theorem. 

Theorem 2. Algorithm 4 terminates. 


Proof. The termination of the Algorithm 4 is based on the reduction of the complexity 
comp (M). As one can see from the comments in the Algorithm 4, after each recursive step 
comp(M') of the new motion polynomial M' strictly decreases. Furthermore, Lines 17-20 
can not happen continually because of /3 < 7 in each motion polynomial. Then in finitely 
many steps we can reduce a and /3 to zero. After this the algorithm will terminate in one 
step using Algorithm 1 . □ 
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3.3. A comprehensive example 

Now we illustrate Algorithm 4 by a comprehensive example where we really enter each 
sub-branch once. We wish to factor the motion polynomial M = P + eD where 

P= (t 2 + 2t + 2)(f 2 + l) 2 , 

D = -(i1 2 + 21 + 2)i + (t 5 + t 4 + 2t 3 + t 2 - t - l)j + (f 4 + t 2 - 2t - l)k. 

First iteration: The input to Algorithm 4 is MW = P W + eQ^ where P^ 1 ) = P and 
lA 1 ) = D from (7). We compute 

Pi = GRPF(P (1) ) = P (1) , T= 1, comp(M (1) ) = (2, 6 , 6 ). 

Thus, we have to use the branch in Lines 32-41 of Algorithm 4: 

hi = —1 — i, h r = — 1 + i, 

Pi = (t 2 + 1) 2 (Z — i + 1), P r = {t 2 + 1 ) 2 (Z + i + 1), ( 8 ) 

Di = jZ 4 + 2 jZ 2 -(i+j + k)Z-l-i-j, = jt 4 + 2 jf 2 - (i+j + k)t + 1 - i - j. (9) 

Note on computation: 

• We compute one quaternion root hi of Pi by (3). We then have Pi = (t — hi)(t — hi) 
and use polynomial division to find Q and R with D = Q(t — h{)(t — hi) + R. The 
dual quaternion h r is the zero of the linear remainder polynomial R. 

• The polynomials Pi and P r are also computed by polynomial division from 

=P r (t-h r ) and PW= T\(t-~hi). 


A similar computation yields Di and D r . 

The updated values of Q, Li and L r are Q = 1, P; = [Zi], L r = [] where Zi = t + 1 + i. 


Second iteration: The input to Algorithm 4 is M ^ = P^ + eQ ^ where P*- 2 -* = Pp 
p( 2 ) = Di are taken from ( 8 ) and (9). We compute 

Pi = GRPF(P (2) ) = (t 2 + l) 2 , T = t- i + 1, comp(M (2) ) = (0,4,5). 


Because of gcd(Pi, ZA 2 )£)( 2 )) = gcd(Pl,TT) = 1 and T ^ 1, we have to use the branch 
in Lines 17-20 of Algorithm 4. Using (3) and polynomial division, we find 


Pi — t 2 + 2 i + 2 , h — —1 + i — ej — ^ek, 


25 


25 1 


M' — t — ^e(7j + k)t 3 + (2 + i|je + i|ek)t 2 — ^e(3j + 4k)Z + 1 — e(i + |tj — §|k). 


( 10 ) 

The updated values of Q 1 Li, and L r are Q = 1, Li = [Zi], L r = [t — h] where = t + h 
and h is as in ( 10 ). 


Third iteration: The input to Algorithm 4 is = P^ + eQ*- 3 -* where M^ = M' is 

taken from (10). We compute 

Pi = GRPF (P (3) ) = (Z 2 + l) 2 , T = 1, comp(M( 3) ) = (0,4,4). 
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Because of gcd(Pi, £>( 2 )p)( 2 )) = 1 and T = 1, we have to use the branch in Lines 9-15 of 
Algorithm 4. Similar to the first iteration we compute 


P 1 =t 2 + 1, P' = t 2 + 1, hi = fi + fj - fk, h r = - i, 
D’ = - i k )Z 3 + (i - ii + - i|k)f 2 


( 11 ) 


_/16_8_ 

v. 35 


35 35 J 


#§j + M k )t 


24 _ 67* 

35 35 J 


_ 6£; _ ALi _ J3_lc 

35 1 175J 175“" 


175 


175 J 


The updated values of Q, Li, and L r are Q = t 2 + 1, Li = [Zi, Z 2 ], L r = [ 72 , 73 ] where 

r 2 =t + i. 


h = t - f i - f j 


A 

7 K > 


Fourth iteration: The input to Algorithm 4 is A/l 4 ) = P( 4 ) + eD ( 4 ) where P^ = 
P'(t — hi)(t — h r ) and = D' are taken from (11). We compute 

Ri = GRPF (P^) = f 2 +l, T = t 2 -(|i-fj+fk)t+f+fj+fk, comp(MW) = (0,2,4). 


Because of gcd(Pi, D^D^)) = 1 and gcd(Pi, TT) = t 2 + 1, we have to use the branch in 
Lines 23-29 of Algorithm 4. Similar to the first iteration we compute 


Pi = t 2 + 1 , P' = t 2 - (f i - f j + f k)i + f + f j + f k, 


hi = - 


158: 


218 : 

483 J 

+31/4 


4011, 
483 
56 


_/ 28 _ 

V75 


J4 j 
575 


483 ‘ 

14: _ Akip _i_ I'd._ 

25 J 25 ' ^ 4" ^ 6g 575 - 

428 : J_ 2096 _ 2308 


h r = — k. 


196 : 

345 J 


7 k )t 2 


345 J 


1725 


1725 


2069 : 
1725 


r 345 J 
613 : , 5531 
1725 J " r 575 


( 12 ) 


The updated values of Q, Li, and L r are Q = (t 2 + l) 2 , Li = [Zi, Z 2 , Z 3 ], L r 
where 


/ j- 1 158 ; _i_ 218 2 _i_ 401 l. „ j. _i_ i, 

I 3 ~ t + 483 1 + 483 J + 483^’ Tl ~ * ■ 


483 


483 


[r- 3 ,r 2 ,ri] 


Fifth iteration: The input to Algorithm 4 is = P^ + eD^ where P^ 5 -* = (t — 
hi)P'{t — h r ) and = D' are taken from (12). Because of Pi = 1, we have to use 
Line 2 of Algorithm 4 and can compute a factorization of M ^ by means of Algorithm 1. 
Because of M^il/P 5 ) = ( t 2 + l) 4 , the factorization is unique. We find M ^ = / 1 / 2 / 3/4 
where 


f _ + _ 158: _ 218: _ 401 1 ,_29- _ 37:, , 2 . 

Jl — b 483 483 J 483 280 e 56 Jt 4 _ 5 tlv ! 

/2 = Z+fi + fj-fk+§d-^je+§ek, 

/ 3 = t - i - fek, / 4 = t - k - fei ■ 


Id- 


Algorithm 4 terminates and the polynomial QM is the product of the ten linear factors 

h, h, h, fi, h, h, U, n, 7 2 , r 3 . 

54. Degree bound of Q 

An upper bound on the degree of Q as returned by Algorithm 4 can be read from the 
following theorem. This degree bound is already know to be optimal. It is attained by 
certain planar motions [1]. 

Theorem 3. The degree of Q as returned by Algorithm 4 is less or equal to the degree of 
the GRPF of the primal part of M. 
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Proof. The proof follows from a careful inspection of Algorithm 4. The increase of the 
degree of Q happen either in Lines 12-13 or lines 26-27. Furthermore, the increase of the 
degree of Q and the decrease of the degree of the GRPF are equal at these places. □ 


We illustrate Theorem 3 by one further example. One achieves the upper bound of 
Theorem 3, the other does not. 

Example 4. The first example is the general Darboux motion considered in [6]. Let M = 
fP — i rjeP £ BH[t] with 


Z = t + 1 ’ *?=2 t_ 4’ P = t ~ h and h = g 1 - gj + g k - 

As seen in [6], this give us the factorization M = Q1Q2Q3, where 

^ 7. 4. 4, 5 . 43 . 97 , 

U 1 = t -1-1 H—k-ei H-ei-ek. 

^ 9 9 9 4 64 J 64 

^ 7. 4. 4, 

Q 2 -t+- 1 +- J --k, 

^ 7. 4. 4, 5 . 43 . 97 , 

Q3 = i -9 1+ 9 J -9 k ^4 el -64 eJ+ 64 ek ' 


Here, no multiplication with a real polynomial is necessary. 

Example 5. The second example is the vertical Darboux motion which was avoided in [6]. 
Let M = £P — irjeP £ 1DHI[£] with 


f = t 2 + 1, 



P = t- i. 


As seen in [ 6 ], no factorization of the shape M = Q1Q2Q3 with linear motion polynomi¬ 
als Q 1 , Q 2 , Q 3 exists. However, we can find a factorization by multiplying with a real 
polynomial whose degree equals the degree of £, the greatest real polynomial factor of the 
primal part of M. We have (t 2 + 1 )M = Q 7 QIQ 5 Q 4 , where 

3 5 3 

Q 7 = t- j --ek, Qe = f+j - -ei + -ek, Q 5 = t- j, Q± = P = t - i. 


3.5. Factorizations in planar motion groups 

Algorithm 4 can produce non-planar factorizations for planar motion polynomials. 
This is an interesting feature but may not always be desirable. If one wishes to find a 
factorization (t — h\) ■ ■ ■ , (t — h n ) of a motion polynomial in a planar motion group, say 
(1, i, ej , ek), with rotation quaternions hi, ..., h n in that group, we have to pick suitable 
left and right factors hi and h r in Algorithm 4. 

Note that for a planar motion in the subgroup (1, i, ej, ek), the primal part and the 
dual part of a motion have a certain commutativity property. If P is a polynomial with 
coefficients in (l,i) and D is a polynomial with coefficients in (ej,ek), then PD = DP , 
e.g., (f — i)ej = ej(t + i) or [t — i)ek = ek (t + i). This allows to transform right factors into 
left factors and vice versa. Moreover, from Equation 3 it follows that there are exactly 
two roots of a real irreducible quadratic polynomial Q in the planar motion subgroup. We 
have, for example, Q = t 2 + 1 = (t — i)(t + i) = (t + i)(i — i). Thus, whenever we compute 
a quaternion root of a quadratic irreducible polynomial in Algorithm 4, we should select 
a solution in the planar motion group and whenever we transfer a left factor hi to a right 
factor h r we should do it in such a way that h r = hi. This ensures that Algorithm 4 really 
returns a planar factorization. 
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